Simple Model for Attraction between Like-Charged Polyions 
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We present a simple model for the possible mechanism of appearance of attraction between like 
charged polyions inside a polyelectrolyte solution. The attraction is found to be short ranged, and 
exists only in presence of multivalent counterions. The attraction is produced by the correlations 
in the condensed layers of counterions surrounding each polyion, and appears only if the number of 
condensed counterions exceeds the threshold, n > Z/2a, where a is the valence of counterions and 
Z is the polyion charge. 



Polyelectrolyte solutions and charged colloidal suspen- 
sions present an outstanding challenge to modern statis- 
tical mechanics. One of the reasons for the great dif- 
ficulty in achieving an understanding of these complex 
systems is due to the intricate role played by the counte- 
rions (microions) which counterbalance the much bigger 
charge of the polyions (macroions). Since the number 
density of counterions is so much greater than that of 
polyions it is the counterions that dominate the ther- 
modynamic properties of polyelectrolyte solutions at low 
densities. The heterogeneity combined with the long- 
range Coulomb force makes polyelectrolytes almost im- 
possible to study by the traditional methods of liquid 
state theory. It is, however, exactly this complexity that 
is responsible for the richness of the behaviors encoun- 
tered in polyelectrolyte solutions and charged suspen- 
sions. 

One of the most fascinating results of the subtle inter- 
play of various interactions is the appearance of attrac- 
tion between two like-charged macromolecules inside a 
solution or a suspension. This attraction is purely elec- 
trostatic and is not a result of some additional short- 
range van der Waals force. The attraction has been ob- 
served in simulations of strongly asymmetric electrolytes 
jjj, as well as a number of experiments The ex- 
act mechanism responsible for this unusual phenomena is 
still not understood, although some theories attempting 
to explain its basis, have been recently advanced |5]]4|]. At 
this time, we believe, there exists a great need to explore 



various simple models which might help shed additional 
light on the origin of this attraction. 

Our discussion will be restricted to rodlike polyelec- 
trolytes, a good example of which is an aqueous solu- 
tion of DNA segments. Let us consider the simplest 
model of such a polyelectrolyte solution. The polyions 
shall be represented by rigid cylinders of net charge — Zq 
distributed uniformly, with separation b, along the ma- 
jor axis of a macromolecule. The counterions shall be 
treated as small rigid spheres of charge aq located at 
the center. The counterions can be monovalent, diva- 
lent, or trivalent, for simplicity however, we shall restrict 
our solution to contain only one of the above types of 
counterions. An appropriate number of coions is also 
present in the solution to keep an overall charge neutral- 
ity. The solvent shall be modeled by a uniform medium 
of a dielectric constant D. It has been argued by Man- 
ning that in the limit of very large Z and infinite dilu- 
tion, a certain number of counterions will condense onto 
the polyions, thus renormalizing their effective charge. 
From a simple phcnomcnological argument Manning de- 
termined the number of condensed counterions to be 
n c = (1— l/a£)Z/a for £ > 1/a and n c — for £ < 1/a, 
where £ = q 2 /Dk B Tb ||. Note that this result subtly de- 
pends on the order of the limits to be taken, first Z — > oo 
and then the infinite dilution. If the limits are inter- 
changed no condensation will appear. Recently we have 
extended the Manning theory to finite concentrations and 
finite polyion sizes. In this case it is possible to show that 
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the counterion association still persists, however, instead 
of a fixed number of condensed counterions associated 
with each polyion, we now find a distribution of clusters, 
each composed of one polyion and 1 < m < Z/a associ- 
ated counterions. The distribution of cluster sizes is well 
localized, and in the limit of large polyion charge and 
infinite dilution approaches a delta function centered on 
the value proposed by Manning, n c 0. 

Can the condensed layer of counterions be responsi- 
ble for the observed attraction between two like charged 
polyions? To answer this question we propose the follow- 
ing simple model. Consider two parallel, rodlike polyions, 
with Z monomers, inside the polyelectrolyte solution. 
The separation between two macromolecules is d. If the 
attraction is produced by some sort of charge-correlation 
mechanism, we expect that it should be short-ranged. 
We shall, therefore, restrict our attention to distances 
such that d < where £d is the Debye screening 
length. As was mentioned earlier, the strong electro- 
static attraction between the polyions and the counte- 
rions favors the formation of clusters composed of one 
polyion and some number of associated counterions. For 
the purpose of this exposition, we shall neglect the poly- 
dispersity of cluster sizes and assume that both polyions 
have n < Z/a condensed counterions. It is important to 
remember that we are concerned with the interaction be- 
tween the two polyions inside a polyelectrolyte solution. 
As was stressed before, an isolated polyion can confine 
counterions only if it is extremely long, while inside a 
solution the cluster formation can take place with the 
polyions of any size . 

The associated counterions are free to move along the 
length of the polyions. We define the occupation vari- 
ables cry, with i — 1,2, Z and j = 1, 2, in such a way 
that cry = 1 if a counterion is attached at i'th monomer 
of the j'th polyion and er^ = otherwise. Since the num- 
ber of condensed counterions is fixed by thermodynamics 
jjj , the values of occupation variables obey the constraint 

52i=i a H = Tlii=i a i2 = n - We shall assume that the only 
effect of the counterion association is a local renormal- 
ization of the monomer charge. The Hamiltonian for this 
model takes a particularly simple form, 



sates for its abstract nature. Evidently, if we can un- 
derstand under what conditions the attraction can arise 
between two macromolecules in this idealization, in the 
future it might pave the way to a more complete, physi- 
cally realistic model. 

With this disclaimer in mind, we proceed to rewrite 
the Hamiltonian as (3TC = £H, where j3 = l/fc^T and £ 
is the Manning parameter defined earlier. The adimen- 
sional reduced Hamiltonian is 



(1 - a(Xij)(l - aai>j>) 
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Using a transformation of occupation variables defined 
by (Tj . = 1 — Uij it is easy to see that the Hamiltonian 
exhibits the following symmetry 

H(Z, a, {a}) = (a- 1) 2 H{Z, a', {a'}), 

where ml = Z — n and a 1 = 1 + l/(ot — 1). 
The partition function is 



(4) 



where the prime indicates that the occupation numbers 
are subjected to the constraints of the number of con- 
densed counterion conservation. The symmetry relation 
(Q) leads to invariance of the partition function, and any 
thermodynamic quantity derived from it, Q(Z, n, £, a) — 
Q (Z, Z - n, [a - 1] 2 £, 1 + 1/ [a - 1] ) . This property is quite 
useful when performing the calculations, since it reduces 
the ranges of parameters needed to study. 

It is convenient to rewrite the partition function in 
terms of the variables yi — exp(— = 1,2, Z — 1, 
associated with the intrapolyion interactions and Z{ — 
exp(— + i 2 ), i = 0, 1,...,Z — 1, related to the in- 
terpolyion interactions. The partition function, for given 
values of Z and n, may be expressed as 
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M — 1 j,3 =1 



q 2 {l - auj 3 ){l - cteyy) 
r(i,j;i',f) 



(1) 



to + {i',f), 



where the sum is restricted 
r(i,j;i',j') — byj\i — i'\ 2 + (1 — Sjj>)x 2 is the distance 
between the monomers located at and (i',j'), Sjj> 

is the Kronecker delta, and x = d/b. Clearly, the above 
model is a great over-simplification of physical reality. 
The molecular nature of the solvent is not taken into ac- 
count. The counterions are assumed to be confined to the 
surface of the polyions, while the polyions themselves are 
treated as completely rigid. Nevertheless, we believe that 
the simplicity of this phenomenological model, which al- 
lows us to perform exact analytic calculations, compen- 
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where the sum runs over all the allowed configurations 
of counterions. The exponents and Vik are quadratic 
polynomials in a with integer coefficients. For not too 
big values of Z, it is possible to generate the whole set of 
integers in the expression above, thus obtaining the par- 
tition function exactly. The force between the polyions 
may now be calculated through 



F = 



1 d\nQ 

bfJ dx ■ 



(6) 
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FIG. 1. Force versus distance between polyions for Z = 20, 
a — 2, £ = 2.283 (corresponding to polymethacrylate) and 
n = 5, . . . , 10 (from top to bottom) in the Monte Carlo simu- 
lation. The lines are only guides to the eye. 



We find that at short distances and for a > 1 the force 
between the polyions may become attractive (negative). 
The exact results are in full agreement with the Monte 
Carlo simulations which can also be extended to much 
larger values of Z, Fig. |l| For a — 1 the force is always 
repulsive, which is in full agreement with the experimen- 
tal evidence on absence of attraction if only monovalent 
counterions are present (2). When the polyion charge 
is completely neutralized, n = Z/a, the force becomes 
purely attractive, as can be seen in Figs. |l| and |2[ Fur- 
thermore, the exact solution and the simulations indicate 
that a critical number of condensed counterions is neces- 
sary for attraction to appear, see Fig. |[ 
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FIG. 2. Distance below which the polyion interaction be- 
comes attractive versus the number of condensed counterion 
in the Monte Carlo simulation. Averages are over 100 sets 
and Z ranges from 20 to 200. 
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FIG. 3. Log-log plot near n — Z/2a for the a — 2 data in 
fig. g. The fit exponent is 0.32. 

To simulate the model we use a standard Monte Carlo 
with particle-hole exchange, not restricted to nearest- 
neighbor pairs. This non local diffusion comes from the 
fact that the counterions may associate/dissociate at any 
point along the polyion. After thermalizing, both the en- 
ergy and the force are measured, results being time and 
sample averaged. In Fig.|^ it is clear that the relevant 
parameter is the fraction of neutralized charge, an/Z, 
the curves for several values of Z (here from 20 to 200) 
collapsing on an universal function, deviations occurring 
only for very small Z. 

As expected, the attractive interaction is short-ranged 
and at larger distances the force becomes repulsive. To 
further explore this point we consider the limit of small 
x. In this case the variable z$ vanishes, while for finite 
£ and Z all the other variables yj and Zj remain posi- 
tive definite. If we define v = mim, (ui,o), it is possible to 
rewrite the partition function as Q = ZqW(1 + P), with 
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where we suppose v i0 = v for the first I of the N c configu- 
rations. As x — > 0, the function P vanishes, and we may 
use the approximation Q w Wzq. We then notice that 
v corresponds to the configurations which maximize the 
number of favorable horizontal interpolyion interactions, 



v = Z — 2na 



(9) 



for a > 2. Using the above approximation for the parti- 
tion function we obtain 



F 



1/1 dW £v 
b/3 \W~dx + x~z 



(10) 
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It is not difficult to see that the derivative of W vanishes 
as x — > 0, so that for small x we can Taylor expand the 
first term of Eq. RG. To leading order we find 



/ = bf3F : 



hx, 



(11) 



where h < (as v vanishes) is a function of the remaining 
parameters of the model. We notice that as x — > the 
force vanishes if v — 0, is large attractive if v is negative, 
and is large repulsive if v is positive. In particular, in the 
limit v — *■ 0~ the equilibrium distance is 



1/3 



(2na 



1/3 
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This scaling behavior is also observed in Monte Carlo 
simulations, Fig. ||and|^. 




FIG. 4. Temperature dependence of the force between 
polyions for Z = 20, n = 6, a — 2 and several values of £: 2.28 
(circles), corresponding to polymethacrylate, 4.17 (squares), 
corresponding to DNA and 20 (triangles) for very low temper- 
atures. The inset shows the distance xo at which the attrac- 
tion first appears as a function of £. We see that for divalent 
counterions, £ = 1 serves as a division point for the low and 
high temperatures regimes. Clearly for £ > 1 the attraction 
is driven by the zero temperature mechanism! 

We would like to stress that for a fixed number of con- 
densed counterions the attraction is insensitive to the 
temperature (see figure ^). In fact the force between 
two lines at zero temperature (£ = oo) is almost ex- 
actly the same as at finite temperature, as long as the 
average electrostatic energy of interaction between two 
condensed counterions is greater than the thermal en- 
ergy, a 2 q 2 /Dd > ksT, where d is the average separation 
between the condensed counterions. The mechanism for 
attraction between two like-charged rods is now clear || . 
At zero temperature the counterions on the two rods will 
take on staggered configuration, i.e. if the site of the 
first rod is occupied by a counterion the parallel site 
of the second rod will stay vacant. If the number of 



condensed counterions is above the threshold n = Z/2a 
the favorable counterion-monomer interactions will out- 
number the unfavorable monomer-monomer interactions, 
thus producing a net attraction at short distances. Us- 
ing the threshold value of n we find that the attraction is 
dominated by the zero temperature correlations as long 
as £ > 2/a. On the other hand applying the Manning 
condensation criterion to the threshold number we see 
that the attraction is possible only if £ > 2/a. We thus 
come to an important conclusion: If the attraction exists 
it is produced by the zero temperature correlations. 

In view of the current interest in the possible mecha- 
nisms responsible for the attraction between like charged 
objects it is worthwhile to make some further comments. 
A most common approach used to study this difficult 
problem relies on the field theoretic methodology similar 
to the one developed in Quantum Field Theory to study 
Casimir forces. The partition function is mapped onto 
an effective field theory which is then studied using a 
loop expansion. Due to the underlying difficulty of this 
process the expansion is usually terminated at the first 
loop level which is equivalent to the so called Gaussian 
approximation. The attraction, in this approach, arises 
as a result of the correlations in the Gaussian fluctua- 
tions P). Clearly from the above discussion this is not 
the mechanism responsible for the attraction between the 
two charged rods which was found by us. The problem 
with the field theoretic approaches used up to now is that 
the ground state energy was not properly treated. 

Instead of allowing for a staggered configuration, which 
we found to be responsible for the attraction observed, 
the field theoretic approaches neglect the discrete nature 
of charges and uniformly smear the condensed counteri- 
ons over the polyion. Clearly at zero temperature this 
can only result in a repulsion! The attraction, then, ap- 
pears only as a finite temperature correction, produced 
by the correlations in the Gaussian fluctuations of the 
counterion charge densities on the two polyions, the ef- 
fect which is much weaker than the one found by us to 
be responsible for the attraction. In view of our results 
we must conclude that the currently used field theoretic 
approaches Q must be reexamined. 

We have presented a simple model for the possible 
mechanism of appearance of attraction between the like 
charged polyions inside a polyelectrolyte solution. The 
attraction is found to be short ranged, and is possible 
only in a presence of multivalent counterions. The at- 
traction is produced by the correlations in the condensed 
layers of counterions surrounding each polyion and ex- 
tends uniformly from zero temperature. The attraction 
appears only if the number of condensed counterions ex- 
ceeds the threshold, n > Z/2a. Using the counterion 
condensation theory to estimate the number of associ- 
ated counterions, we see that the attraction is possible 
only for polyelectrolytes with £ > 2/a. This is the fun- 
damental result which has, evidently, gone unnoticed in 
the previous studies of this interesting phenomena. 
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